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Abstract 

This paper discusses the design and performance of the time measurement technique 
and of the synchronization systems of the CMS hadron calorimeter. Time measure- 
ment performance results are presented from test beam data taken in the years 2004 
and 2006. For hadronic showers of energy greater than 100 GeV, the timing resolution 
is measured to be about 1.2 ns. Time synchronization and out-of-time background re- 
jection results are presented from the Cosmic Run At Four Tesla and LHC beam runs 
taken in the Autumn of 2008. The inter-channel synchronization is measured to be 
within ±2 ns. 
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1 Introduction 

The primary goal of the Compact Muon Solenoid (CMS) experiment [H is to explore parti- 
cle physics at the TeV energy scale exploiting the proton-proton collisions delivered by the 
Large Hadron Collider (LHC) 0. The central feature of the CMS apparatus is a supercon- 
ducting solenoid, of 6 m internal diameter, providing a field of 3.8 T. Within the field volume 
are the silicon pixel and strip tracker, the crystal electromagnetic calorimeter (ECAL) and the 
brass/ scintillator hadron calorimeter (HCAL). Muons are measured in gas-ionization detectors 
embedded in the steel return yoke. In addition to the barrel and endcap detectors, CMS has 
extensive forward calorimetry. 

The primary purpose of the HCAL is the measurement of hadronic energy from collisions in 
CMS. In addition to the energy measurement, the HCAL is also able to perform a precise time 
measurement for each energy deposit. Precise time measurements are valuable for excluding 
calorimeter noise and energy deposits from beam halo and cosmic ray muons. Time informa- 
tion can also be valuable for identifying new physics signals such as long-lived particle decays 
and slow high-mass charged particles HO- 

This paper is organized as follows. Section [2] reviews the pertinent details of the HCAL con- 
struction and their fundamental impact on timing resolution. Section [3] discusses the method 
used to extract a time value from the digitized HCAL signal. In Section |4j the validation of 
the method is presented based on measurements in test beam and initial beam operations of 
the LHC in September 2008. Section [5] details the performance of timing filters in the suppres- 
sion of non-collision-based backgrounds and the effect these timing filters have on simulated 
physics events. 

2 CMS Hadron Calorimeter Description 

A detailed description of the HCAL can be found elsewhere d. Briefly, the HCAL consists 
of a set of sampling calorimeters. The barrel [4j and endcap [5I calorimeters utilize alternat- 
ing layers of brass as absorber and plastic scintillator as active material. The scintillation light 
is converted by wavelength-shifting fibers embedded in the scintillator and channeled to hy- 
brid photodiode detectors via clear fibers. The outer calorimeter [6] utilizes the CMS magnet 
coil/ cryostat and the steel of the magnet return yoke as its absorber, and uses the same active 
material and readout system as the barrel and endcap calorimeters. The forward calorimeter is 
based on Cherenkov light production in quartz fibers and is not discussed in this paper, due to 
its different signal time structure. 

The HCAL is segmented into individual calorimeter cells along three coordinates, rj, (p, and 
depth. The (p coordinate is the azimuthal angle and tj is the pseudorapidity. The depth is an 
integer coordinate that enumerates the segmentation longitudinally, along the direction from 
the center of the nominal interaction region. The layout of the barrel, endcap and outer calori- 
meter cells is illustrated in Fig. [I] Most cells include several scintillator layers; for example, in 
most of the barrel all 17 scintillator layers are combined into a single depth segment. 

Calorimetric measurements are acquired using the HCAL readout electronics, shown schemat- 
ically in Fig. [2] Each electronics channel collects and processes the signal of a single cell. One 
calorimetric measurement is acquired with each LHC clock tick (25 ns) from each cell in the 
HCAL. This defines a "time sample". The HCAL front-end electronics does not sample the 
signal instantaneously; rather, the electric current collected from the photodetectors is inte- 
grated over each clock period and then sampled. As a consequence, the sample clock is most 
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2 CMS Hadron Calorimeter Description 




Figure 1: A quarter slice of the CMS HCAL detectors. The right end of the beam line is the 
interaction point. FEE denotes the location of the Front End Electronics for the barrel and the 
endcap. In the diagram, the numbers on the top and left refer to segments in 77, and the numbers 
on the right and the bottom refer to scintillator layers. Colors/ shades indicate the combinations 
of layers that form the different depth segments, which are numbered sequentially starting 
at 1, moving outward from the interaction point. The outer calorimeter is assigned depth 4. 
Segmentation along (p is not shown. 
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Figure 2: A schematic view of the HCAL front-end readout electronics. The readout for one 
HCAL cell /channel is shown. Key features are the optical summing of layers, charge integra- 
tion followed by sampling and digitization, and per-channel programmable delay settings. The 
"QIE" [7\ is a custom chip that contains the charge-integrating electronics with an analog- to- 
digital converter (ADC). The Configuration Data input defines the sampling delay settings. 
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commonly termed the "integration clock". 

The integration clock can be delayed with respect to the LHC clock on a per-channel basis by 
programmable settings, referred to as sampling delay settings, that have a resolution of 1 ns. 
The purpose of these settings is to compensate for channel-dependent timing variations at the 
hardware level, and to permit the energy estimate from each LHC crossing to be reconstructed 
consistently for use in the Level-1 trigger for all pulse amplitudes and for all channels. They 
also provide an initial coarse timing calibration. 

2.1 Sources of timing variation and uncertainty 

There are two dominant sources of channel-dependent timing variation: the different time- 
of-flight of particles from the interaction region to each HCAL cell and the different signal 
propagation times through optical fibers of different lengths. Within the barrel, the first effect 
tends to skew reconstructed times later for higher tj. The second effect, because of the location 
of the front-end electronics, tends to skew reconstructed times earlier for higher t], and this is 
the effect that dominates. The combination of these two effects induces a timing dependence 
on rj, with a spread in each half-barrel of about 15 ns. Smaller variations as a function of (p 
are induced by clock distribution differences and other effects. By applying proper sampling 
delay settings, this spread can be reduced substantially. In addition, the reconstruction software 
utilizes a set of per-channel calibration constants to synchronize the mean timing of energy 
measurements to less than 1 ns. 

Fiber lengths also differ within each calorimeter cell along the radial coordinate, but in this case 
there are no means available to compensate for these differences. Since the signals from all the 
scintillator layers comprising one cell are optically summed, and the optical path lengths are 
not equalized across the layers, the resulting signal is smeared in time. For hadrons, which 
exhibit large shower-to-shower fluctuations in longitudinal development, the optical summing 
of layers imposes a limit on the timing resolution, estimated to be approximately 1 ns. For sig- 
nals with uniform energy deposit in each layer (such as those arising from the beam-collimator 
interactions described later), the resolution is not limited by shower-development fluctuations 
and can be significantly smaller than 1 ns. 



3 HCAL Time Reconstruction 

When the Level-1 trigger system [8j identifies an event of potential physics interest, a set of 
10 consecutive time samples per channel containing the triggered bunch-crossing is collected 
and sent to the high-level trigger software. The capability of the HCAL system to reconstruct 
the arrival time of the signal more precisely than the sample clock period derives from the 
spread of the HCAL pulse shape over three to four time-integrated samples (Fig. [3]). The time 
reconstruction software calculates a first order time estimate from a center-of-gravity technique 
using the three samples centered on the peak, 



Weighted peak bin 



(p-l)Ap_ 1 + pA p + (p + l)A pH 
Ap-i + Ap + Ap+i 



xC , (1) 



where A, represents the amplitude of an arbitrary time sample i, and p is the value of i such 
that Ap is maximum over the set of samples. In the case of multiple samples with the same 
amplitude, the earliest one is picked. The constant C is an amplitude-independent normaliza- 
tion constant that rescales the first order estimate to a range from zero to one. The weighted 
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4 Detector Timing/Synchronization Commissioning 



peak bin is then used to determine a second order correction (Fig. |4]| that compensates for the 
asymmetry of the pulse shape, yielding the phase of the signal within the peak time sample. 
An additional additive correction determined in test bench measurements 0, 



Asiew = (3.59 ns) log 



10 



ITeV 
E(TeV) 



(2) 



compensates for an electronics effect that delays the measured time for signals with pulse en- 
ergy E less than 1 TeV. The maximum value of A s i ew is taken to be 10 ns. 

This algorithm yields accurate time measurements except when consecutive bunch crossings 
yield energy measurements of similar amplitude within the same HCAL cell. Such events can 
be identified by their anomalous pulse shape, and the time measurement can be marked as 
having poor resolution. 



Pulse Shape 

Integrated Pulse 




-10 1 2 3 4 

Sample Clock Time (bins) 




0.6 0.7 0.8 0.9 1 

Weighted Peak Bin 



Figure 3: The HCAL pulse shape and its re- 
lation to the integrated samples. Time sam- 
ple is defined by the trigger for a nominally 
synchronized event. 



Figure 4: The correction function applied 
in signal reconstruction software to compen- 
sate for pulse shape asymmetry. The straight 
line represents a trivial correction function 
for a perfectly symmetric pulse. 



4 Detector Timing/Synchronization Commissioning 

This section describes the performance characterization and commissioning results for HCAL 
timing and synchronization. Beam tests at the CERN SPS H2 area (hereafter referred to as "test 
beam") produced key results [10] that define the timing performance of the HCAL barrel and 
endcap; these results are presented first. The results from the LHC 2008 beam commissioning 
data for the barrel, endcap and outer calorimeters are then described; these results validate and 
extend the initial synchronization calibration constants determined from beam tests. 

4.1 H2 beam test data 

For test beams in 2004 and 2006, a section of the HCAL barrel was mounted on a table that was 
adjustable in rj and cp coordinates. The unit under test was exposed to pion beams with energies 
in the 3-300 GeV range. The delivered beam was asynchronous with the 40 MHz integration 
clock, so all relative sampling phases could be investigated. From these tests the following 
important timing benchmarks for the HCAL were accurately measured. 



4.1 H2 beam test data 
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Figure 5: Average pulse shape reconstructed 
for the HCAL from H2 beam test data using 
300 GeV pions. 



Figure 6: Resolution performance of the 
HCAL barrel from 150 GeV pion test beam 
data. The non-Gaussian tail is attributed 
to beam contamination from electrons and 
photons that shower in the first layers of the 
HCAL and thus reconstruct with later signal 
times. 



The HCAL pulse shape was successfully reconstructed from data. Since the pulse is integrated 
before it is digitized, this required a numerical differentiation technique that utilized the phase 
between particle arrival and the integration clock, as measured by a time-to-digital converter 
(TDC) running at 32 times faster than the sample clock rate. This allowed the separation of 
the data sample into subsamples S„, with n =1.32 ordered by TDC phase. These subsamples 
represent the integration of the pulse P(t) shifted successively by the TDC bit resolution (At = 
0.78 ns). Then, 



rT+At \/J\ 

(Im) Sn+1 ~ (Im) Sn * I P(t)dt~ Pavg(T)A£ , P avg (r) c± -A_> 



T 



(3) 



where {I m )s„ is the average integrated amplitude for time bin m of data subsample S n , x = 
mAt is the phase within that time bin, and P a vg(i") is the average value of P(t) over At at 
t = t. Equation [3] is valid to reconstruct the first 25 ns of the pulse shape, for which m = 0. 
For each remaining 25 ns interval, bin m is incremented and a second integral term appears 
in the equation that corresponds to the Ai-sized portion of the pulse "leaving" the time bin. 
In this manner, the whole pulse was reconstructed to approximately 0.8 ns resolution, as is 
shown in Fig. [5] The functional form matched to these results is used in the simulation of the 
CMS detector response, and also to determine the second-order time reconstruction correction 
function (Fig. [4]). 

To measure both the channel-to-channel variations and the hadronic timing resolution limit of 
the HCAL, a 150 GeV pion beam was directed at every cell in the unit under test. From these 
data the variation of the timing of signals as a function of rj between HCAL cells in the barrel 
was mapped. As the test setup mimicked the geometry and fiber lengths of the final barrel 



detector, a set of final sampling delay settings could be derived. Section 4.2 discusses the per- 
formance of these settings with LHC beam data. The HCAL timing resolution was determined 
by removing the channel-to-channel variations and combining the data. A Gaussian function 
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was then fitted to the distribution of the leading energy deposit from each event, producing 
a best-fit width of a = 1.2 ns (Fig. [6]!. This result is consistent with the expected effect of the 



uncorrectable sources of timing spread described in Sec. 2.1 averaged over many hadronic 



shower events. The measured resolution is consistent across the full set of channels studied. 

The variation of HCAL timing resolution as a function of energy was also measured with pion 
beam data. In this measurement the HCAL unit under test was paired with its correspond- 
ing section of the ECAL in front, as in the CMS detector. Events in two categories were con- 
sidered: those with no interaction in the ECAL (<lGeV energy deposited) and those with a 
significant interaction (>5GeV energy deposited). The results indicate that the presence of 
the ECAL does not substantially alter the timing resolution of the HCAL as a function of en- 
ergy (Fig. [7](a)). CMS software was subsequently adjusted by smearing the times of simulated 
HCAL energy deposits so that simulated data of the test beam configuration exhibited the same 
energy-dependent time resolution (Fig. [7](b)). 



"a 

CO 
CD 

CL 
CO 



10 



-10 



[CMS 



TB06 ECAL/HCAL Energy Scan 

[ ' ^ Non-interacting in ECAL (95% 

| | Non-interacting in ECAL (68% 

Interacting in ECAL (95%) 

Interacting in ECAL (68%) 



50 



(a) 



100 150 200 250 



300 
"cell 



350 



E„„ (GeV) 



C/3 

"O 

CO 
CD 

Q_ 

CO 



10 







I, 


CMS 




Simulation, Pion Gun 

Smeared Timing 
I I 95% 

68% 

Time Filter 






















































































(b) 








^ , , , 


, , , , 


, , , , 











50 



100 150 200 250 



300 
"cell 1 



350 



E„,, (GeV) 



Figure 7: (a) Timing resolution as a function of energy measured during test beam runs in 2006, 
showing the consistency of time reconstruction for particles that begin showering in the ECAL 
(lines) and those that do not (areas), (b) Timing resolution as a function of energy of recon- 
structed deposits generated from CMS simulations. The lines marked "Time Filter" indicate 
the boundaries of a timing window discussed later and used to accept or reject deposits for 
distinguishing signal from background. 



4.2 LHC beam data 

Data from LHC beam commissioning were used to validate the barrel sampling delay settings 
derived from test beam data, and to derive settings to compensate for the rj dependence of 
timing in the endcap and outer calorimeters. Data from LHC beam operations were recorded 
for the first time on September 10, 2008. In preparation for the arrival of LHC beam at CMS, 
the //-dependent sampling delays measured from test beam were loaded into the barrel front- 
end electronics. These delay values had been calculated to synchronize channels for data from 
collisions. Delay values for endcap and outer sectors were set to zero, since no timing measure- 
ments for these sectors were available at the time. 

The LHC beam delivered to CMS consisted of a single 450 GeV proton bunch circulating in 
the ring, first in one direction and then in the other. The LHC beam was also steered into 
collimators located 150 m upstream in either direction of the detector interaction point. Such 
events are referred to as "beam splash" events. Because of the large amount of material in 
the collimator and along the path of the secondaries, only muons were able to penetrate the 
entire CMS detector and form the signal in the calorimeter system. The millions of muons 
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produced from the dumping of the beam passed through the entire detector, depositing large 
signals in every cell of the HCAL, with equivalent energy ranging between hundreds of GeV 
to a few TeV. A schematic representation of the geometry of the "beam splash" setup is shown 
in Fig. § 
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Figure 8: A schematic view of the geometry of "beam splash" events. This geometry was 
used to predict the timing of energy deposits and thereby highlight channels in the HCAL 
that required synchronization. The times Tjp and T co \ refer to the times-of-flight of relativistic 
particles from the interaction point (IP) and collimator, respectively. The timing numbers give 
the specific values of Jn> to the front faces of the barrel and endcap. 

The signal from "beam splash" muons has a well-defined time with respect to the arrival of 
the LHC bunch and to the LHC clock. The time-of-flight of the muons from the collimators 
to the geometrical centers of each cell was calculated and used to predict the expected timing 
versus rj for all HCAL cells. The data were averaged over all (p values at each rj value to 
obtain an HCAL rj timing profile. The resulting distribution of the differences between the 
measured and predicted times is shown in Fig. |9ja). The timing measurements for the barrel 
region are distributed around zero within ±1 ns, demonstrating the correctness of the delay 
settings calculated from beam tests. The timing for the endcap before any phase corrections 
exhibits an average time difference that is 10 ns later than the barrel. The x-shaped pattern and 
the structure at \rj\ > 1.4 in the outer calorimeter are understood to originate from the pattern 
of optical fibers in the scintillator trays and the positioning of readout electronics. 

From these data, //-dependent sampling delays were derived for the endcap and outer calorime- 
ters. The delays were loaded into the front-end electronics, and new "beam splash" data were 
taken from interactions on the +z collimator on September 18, 2008. The same analysis used to 
derive the original delays was repeated on these data; the results are shown in Fig. |9pb). Some 
outer calorimeter data points are omitted from the plot due to a technical issue at the time of 
data-taking that has since been resolved. 

Figure [9pb) shows that, as a result of the sampling delay derivation from "beam splash" data, 
barrel and endcap channel synchronization along rj is verified to within ±2 ns. Although 
"beam splash" event data from both directions were used to derive the settings for endcap 
and outer, only event data from the +z direction were available to verify these settings. For this 
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Figure 9: Difference between measured and predicted time as a function of tj from LHC "beam 
splash" events, (a) Per-channel delay settings that compensate for timing variations along rj 
were loaded into the front-end electronics for the barrel (depths 1-2, \tj\ < 1.4) but not for 
the endcap (depths 1-3, \rj\ > 1.4) or outer (depth 4) calorimeters. Events from both +z and 
— z directions are included, (b) Compensating sample delay settings were loaded for all three 
calorimeters, barrel, endcap and outer. Only events from +z direction were available and are 
included. 

reason Fig. |9^b) exhibits systematic deviations, particularly near the barrel-endcap boundary. 
For cells that are significantly slanted relative to the arrival direction of the muons, the "beam 
splash" will illuminate the various layers of the cell at different times. The effect is canceled 
when considering events from both directions. Additional possible systematic effects, such 
as an offset in timing between the positive and negative endcaps, will be studied with future 
"beam splash" data from both directions. 

The results shown in the last two sections indicate that the methods of time reconstruction 
and synchronization are effective both in the H2 test beam environment and with the HCAL 
integrated in CMS. The final offline time corrections, which will be determined with collision 
data when they become available, are expected to provide synchronization with a spread in the 
per-channel mean times of less than 1 ns. 



5 Impact of Timing in Analysis 

Many interesting physics processes that the CMS experiment was designed to study or discover 
contain the signature of an undetected particle. Examples include the Standard Model neutrino 
and the lightest stable particle in several hypothetical supersymmetric particle spectra. Such 
particles would pass undetected through CMS, carrying some of the energy and momentum 
that would otherwise counterbalance the other collision products registered in the calorimeters. 
At hadron colliders, the initial energy and longitudinal momentum of the colliding hadronic 
constituents are not known, but the initial transverse momentum is known to be negligible in 
comparison and can be treated as zero. Therefore, the benchmark used to infer the existence of 
an unobserved particle in any given event is the calorimetric missing transverse energy [11J. 

The missing transverse energy (MET) is reconstructed from calorimeter towers, which are re- 
constructed objects containing the sum of the signals from HCAL and ECAL cells at the same 
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(tj, <p) coordinates. The MET is calculated by taking the magnitude of the vector sum of the 
transverse energy contributions of towers, 



where E, and (0;, (pi) are the energy and angle coordinates of each tower. 

Multiple complications arise in the accurate measurement of MET. The energies of jets or lep- 
tons can be mismeasured due to detector miscalibration or shower fluctuations in the dead 
material of the detector. The MET in a collision event can also be affected by energy deposits 
in HCAL cells from non-collision sources. Without a separate identification of these sources 
the energy will be misassigned to the event of interest, inducing fake MET. Example sources 
include cosmic ray muons, detector noise, and beam halo. 

This section explores the use of HCAL timing as a potential tool to identify sources of fake 
MET. The CMS design ensures that the phase of collision products relative to the LHC clock 
is fixed to high precision. The results presented in Section |4] indicate that the HCAL system is 
capable of being synchronized to approximately 1-2 ns. Since sources of fake MET are generally 
uncorrelated with the collision event one can significantly suppress their contributions to MET 
using their reconstructed times. In the case of beam halo particles, although their arrival is 
correlated to the arrival of proton bunches in the beam, the exact arrival time to a given detector 
cell is generally distinguishable from the arrival time of particles produced in the interaction 
region. 

5.1 Impact on out-of-time background 

The strategy for using timing to reject backgrounds is implemented at the level of reconstructed 
cell measurements, rather than higher-level objects like hadronic jets or MET. It is expected 
that during data taking triggered events will contain collision data in around 1000 HCAL cells, 
over which additional out-of-time energy deposits may be superimposed. These out-of-time 
deposits can be individually removed from the calculation of MET rather than rejecting the 
entire event. The application of such a technique is appropriate for the analysis of prompt 
physics signals, while an analysis focused on long-lived particles might perform the reverse 
selection, preferring energy deposits that are out of time. 

The sources of irreducible timing uncertainty described in Section 2 will limit the tightness 
of the timing requirements that can be applied. The degradation of time resolution for low 
energy particles requires that any time filtering must take this energy dependence into account 
in order to avoid biasing the measurement. The timing filter for individual cell measurements 
was defined according to the resolution measurement shown in Fig. [7| The filter, shown in 
Fig. |7jb), is made wider than the 95% confidence level by a factor of two, to allow for the 
additional uncertainty associated with the spread of primary vertex locations, particularly for 
the endcap calorimeters. The effect of the primary vertex distribution is expected in collision 
data and included in the CMS simulation but not included in Fig. |7jb), which is intended to 
replicate the results of test beam data. The larger envelope is also considered a conservative 
choice for startup when the detector is expected to be less well calibrated, and other sources of 
irreducible uncertainty will not yet have been identified. For hits with energies below 4 GeV, no 
requirement is applied except that the peak amplitude should be in the triggered data sample. 

Using the strategy outlined above, the timing filter was then applied to two datasets: beam 
background data collected during LHC beam commissioning in September 2008 and cosmic ray 
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Figure 10: Time filter performance on MET for HCAL-triggered noise events, events selected 
by the cosmic ray muon trigger, and events from beam operations. The top plots show the 
MET distributions calculated twice for each event, once using all measurements (unfiltered) 
and once using only in- time measurements (time-filtered). The bottom plots show the fraction 
of events that survive as a function of a MET threshold applied to both unfiltered and time- 
filtered measurements. 
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muon and CMS detector noise background data collected in October and November 2008, dur- 
ing the month-long data-taking exercise known as the Cosmic Run At Four Tesla (CRAFT) [12J. 
For each event, the MET was calculated twice for comparison. One calculation of MET used 
all energy deposits (referred to as "unfiltered"), while the second did not include energy de- 
posits with a measured time outside the window (referred to as "time-filtered"). The unfiltered 
measurements include energy depositions in bunch crossings other than the triggered one, as 
a consequence of the triggering and reconstruction algorithms used during these periods. 

The comparison between the unfiltered MET and the time-filtered MET for HCAL-triggered 



events, cosmic ray events, and beam events is shown in Fig. 10 The HCAL-triggered events 
were collected using a Level-1 trigger in which a single calorimetric tower energy exceeds a 
5 GeV threshold. Therefore, this sample includes photo-detector and other instrumental noise 
events |Tl3| . Events with a Level-1 muon trigger were removed from the sample. The cosmic 
ray muon events were collected during CRAFT with a Level-1 trigger in which a single muon 
is identified by the CMS muon stations. The beam events were taken with a single circulating 
450 GeV proton bunch on September 18, 2008. All events from this sample are included in the 
plots regardless of trigger source, which could be beam pickup, coincidence of scintillator or 
forward calorimeter energy deposits, or any of several muon triggers. The beam sample thus 
includes beam halo and beam-gas events as well as cosmic ray muons and accidental overlaps 
with calorimeter noise events. In all cases, the physical MET value would be zero, so a net 
migration of events from high MET to low MET is expected when comparing the unfiltered 
and time-filtered MET distributions. 



Figure 10 also shows the fraction of events that survive a minimum MET criterion as a function 



of the threshold value. This was done for both the time-filtered and unfiltered distributions 



for each background type. Figure 10 therefore, indicates that a reduction by a factor of 5-10 
in the rate of events with fake MET can be achieved with time filtering for all three sources 
of background studied. This rejection factor is consistent with the ratio of the filter window 
width to the width of the four time samples considered in the reconstruction of an HCAL 
measurement. The fake MET production for all three samples is due predominantly to a small 
number of high energy measurements, against which the filter is most effective. 

The filter can be tuned to enhance further the background rejection performance. More sophis- 
ticated algorithms could be employed to filter clusters of cells with mixed in-time and out-of- 
time signals if the members of the cluster share common characteristics that identify them as 
background. For instance, beam halo muons may create a cluster of deposits in consecutive 
cells along tj for the same (p, most of which will be out-of-time. This pattern, once identified, 
could allow the removal of nominally in-time background as well. 

5.2 Impact on physics processes 

To estimate the impact of this technique on signals that have not yet been produced in the CMS 
detector, studies have been performed on two simulated processes, one that contains MET at 
the event generator level and a second that does not. An interesting example of the first case 
is the real MET induced by the hypothetical weakly-interacting supersymmetric (SUSY) |fl4| 
particle called the neutralino. By applying an ill-chosen time filter, one may underestimate 
the MET induced by the neutralino by filtering out calorimetric deposits associated with the 
recoiling jet. A typical process without significant MET at generator level is a hadronic dijet 
event, commonly referred to as a QCD dijet event. In this case one may also infer a (fake) MET 
in an event in which no MET is expected, causing a tail in the high end of the MET distribution. 
Therefore, an ill-chosen time filter could degrade the discovery potential of a neutralino search 
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by reducing the measured MET of the neutralino signal while simultaneously increasing the 
(fake) MET of a major background. 

To assess whether the chosen timing cuts induce a bias on processes with significant MET, the 
timing filter was applied to a sample of simulated events containing neutralinos and generated 
under the assumptions of \/s = lOTeV and no overlapping events (L = 10 30 cm _2 s _1 ). The 
selected SUSY sample was generated using ISASUSY IIT5H and Pythia 6.4 ffT6| at a "low mass" 
point referred to as "LM1" fF7\ in the phase space of one favored model of SUSY known as 
mSUGRA (Mj. The early discovery or exclusion of this point is potentially within reach of 
the LHC experiments given the integrated luminosity expected from the first physics run. The 



left plots in Fig. 11 show that the same energy-dependent timing window used in Section 5.1 



preserves MET distributions for simulated SUSY events with high efficiency. 

The impact of the timing filter for QCD dijet events was also studied. The QCD dijet process 
was generated by Pythia for a range of hard scattering scales 15GeV/c < pj < 3TeV/c with 
same center-of-mass and luminosity assumptions as the SUSY sample. While these events 
have no inherent MET, mismeasurements and dead material will result in the reconstruction 
of some MET in these events, particularly for events with large total transverse energy. The 



results for the QCD dijet sample are shown in the right column of Fig. 11 and indicate that the 
MET distribution is maintained by the filtering process. It is important to note that this filter 
is not intended to remove QCD backgrounds, which are inherently in-time. The point of the 
comparison is to verify that an increase in MET is not inadvertently created by the filtering 
process. 

In order to apply the time filter to analysis of physics processes, actual performance of the filter 
will have to be rederived on collision data using processes such as W — > iv events containing 
jets. The generation of fake MET will be measured using processes such as 7 + jet, where no 
MET is expected and the recoil against the photon provides good control of possible jet energy 
mismeasurement. The efficiency of background rejection can be further studied using dijet 
event samples containing a cosmic ray or beam halo muon as identified by the muon system. 



6 Summary 

This paper has presented the technique used to determine the arrival time of energy deposits 
in the CMS HCAL to less than 25ns, and the use of this technique to suppress out-of-time back- 
grounds. The performance of the technique and of the synchronization system of the HCAL 
have been demonstrated using data from test beam operations as well as "beam splash" data 
from the LHC in September 2008. The ultimate resolution of the technique for hadron show- 
ers was measured to be 1.2 ns using test beam data, and the time spread of HCAL channels 
after hardware alignment in t] was measured to be ±2 ns, with an expected potential offline 
synchronization of less than 1 ns. Using data from CRAFT and beam operations in 2008, the 
use of a timing filter to suppress fake transverse missing energy was studied. This technique 
was effective in reducing the rate of MET produced by cosmic ray muons, detector noise, and 
beam background by a factor of 5-10. At the same time, the integrity of simulated signal MET 
distributions under the same filtering technique was verified. This technique and the timing 
results from HCAL in general are available for use in CMS analyses to suppress backgrounds 
or select particular signals. 
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Figure 11: Time filter performance for SUSY and QCD processes showing the spectrum of fil- 
tered and unfiltered MET (top plots) and the MET survival rates (bottom plots). These indicate 
that the MET is not significantly affected by the operation of the time filter for in-time pro- 
cesses. These simulated samples were generated at = lOTeV, L = 10 30 cm" 2 s _1 , and with 
no pileup. 
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